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We show that the diffractive structure function is perturbatively calculable in the domain 
where the diffractive mass is small but still outside the resonance region. In this domain, 



which can be characterized by A 2 /Q 2 <C 1 — /3 <C yA 2 /Q 2 , the structure function repre- 

sents a new observable, which is highly sensitive to the small-x skewed gluon distribution. 
Our leading order calculation and the estimate of next-to-leading order corrections are 
consistent with available data and demonstrate the potential of more precise data to put 
further constraints on skewing effects. 



Abstract 




1 Introduction 



During the last few years, the concept of skewed parton distributions (see [l|] and refs. 
therein) has attracted much interest. These quantities, which are also known as off- 
diagonal or non-forward parton distribution functions, represent an extension of conven- 
tional parton distributions and allow for a QCD analysis of hard processes where the 
hadronic target is scattered elastically 0. 

In this letter, we focus on the role of the skewed gluon distribution in small-mass in- 
clusive diffraction, i.e. the (3 — > 1 limit of F®(j3, £, Q 2 ). The diffractive structure function 
F 2 D is defined by the process 7*p — > p'X. Our variables are the photon virtuality Q 2 , 
the invariant mass M of the diffractive system X, and the two dimensionless quantities 

P = Q 2 /(Q 2 + M 2 ) and £ = xjp = x B j//3. 

We propose the measurement of the diffractive cross section in the kinematic domain 
Q 2 M 2 3> A 2 (where A ~ Aqcd is a hadronic scale) as a new method for the extraction 
of the skewed gluon distribution at small x. At the same time we point out that, given 
the skewed gluon distribution, the diffractive cross section in this region is perturbatively 
calculable. This corresponds to the calculation of for f3 approaching 1 but still outside 
the resonance region M ~ A. We compare our approach with the limited large-/? data 
available at present and discuss expected next-to-leading oder corrections on the basis 
of the unintegrated gluon distribution. 



2 at very large Q 2 and /3 — > 1 

The following analysis is based on the well-known formulae || for the qq pair production 
by a virtual photon scattering off a hadronic target. At lowest order, the energetic photon 
fluctuates into a qq pair with quark, antiquark momentum fractions z and 1 — z before 
it reaches the target. For a fixed value of z, the typical transverse size of the pair is 
~ 1/e, where e 2 = z(l — z)Q 2 + m 2 with m the mass of the quark ||. If 1/e is sufficiently 
small, the proton target can be characterized by its gluon distribution at the scale e, as 
is known from vector meson electroproduction ||||. 

The relevant formulae for longitudinal and transverse incoming photons producing a 
pair of quarks with electric charge e q read 
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and the kinematical limits of the z integration are determined by 



' (4) 
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The prefactor 1/6 is due to the t integration performed under the assumption of a t 
dependence ~ exp [bt]. 

Equation ([!]) is valid in the leading log(l/x) approximation, where the difference 
between the longitudinal momenta carried by the two exchanged gluons is irrelevant |7| 
and the conventional gluon distribution appears. To take this difference into account, the 
gluon distribution has to be replaced by the skewed gluon distribution H g according to 

G(x, jj 2 ) = H g (x, x, fx 2 ) — »■ H(x,x',/j 2 ), (5) 

where x and x' are the momentum fractions of the proton carried by the first and 
second exchanged gluons. With the main contribution coming from the region x' <C x, 
the resulting correction can be accounted for by introducing the multiplicative factor 
R g = H g (x,x' <C x, fi 2 )/H g (x, x, /i 2 ) into Eq. (1) (cf. the analyses of ||(|). This skewing 
factor has been calculated for x <C 1 in IM to be 



9 v^F T(A + 4) ' dln(l/a;) V ; 



In deriving Eq. the underlying amplitudes were assumed to be purely imaginary. 
This can only be justified in the leading log(l/x) approximation. Therefore, in addition 
to skewing, the effect of the real parts of these amplitudes has to be taken into account. 
Under the assumption of appropriate crossing properties of the amplitude T and of a 
purely power-like behaviour of its imaginary part (i.e. a power-like behaviour of the gluon 
distribution), 

ImT ~ s x (where s = Q 2 /xBj), (7) 

the relation 

ImT V 2 J K J 

can be derived (see, for example, |5|,[TT|). 

Equation ([!]) with the correction factor R 2 (l + C 2 ) evaluated at x = £ and fi 2 = e 2 
forms the basis of our leading order analysis. 

Consider first the simple case when m = and (5 — > 1. In this limit, fi{z) = 1 and 
Zmin = so that Eq. ([]]) provides a parameter-free prediction for the longitudinal cross 
section at asymptotically large Q 2 . This prediction is not affected by the endpoints of 
the z integration since the effective anomalous dimension 7 of the gluon distribution in 
this region, 

xG(x,/x 2 )~(/i 2 ) 7 (7>0), (9) 

suppresses the integrand (cf. |]I2| ). 



3 



The transverse cross section at m = and (3 — > 1 is suppressed with respect to 
the longitudinal one by an explicit factor ~ (1 — (3) contained in fr{z)- However, the z 
integration in the transverse case is divergent even for 7 > 0. Since no free quarks can 
be produced, it is natural to use a cutoff z m i n ~ A 2 /M 2 , which follows from Eq. if 
m is replaced by a hadronization scale ~ A. (Note that for sufficiently large (3 the scale 
^Q 2 always stays in the perturbative domain.) The z integral in 



> Z r , 



then gives a result ~ l/z^^ 7 leading to an overall ratio 



This ratio is small as long as 
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which, even when combined with the condition M 2 ^> A 2 , leaves a window for perturba- 



tive calculability of F£r ■ A particularly simple, conservative expression for this window 



is obtained by setting 7 
leads to 



in Eq. (|TT|) , which on switching to the kinematic variable j3 



A 2 /Q 2 < 1 - p < JA^/Q 



(12) 



In this region, the diffractive structure function is dominated by the perturbatively calcu- 
lable, longitudinal, higher-twist contribution, the phenomenological importance of which 



has been emphasized before [13 



Contributions associated with the production of qqg final states are known to be 
important for F 2 D [14]]. However, they are not relevant in the limit /3 — > 1. This can be 
seen by relating the qqg final state to the diffractive gluon distribution [15[], which is 
known to fall like (1 — j3) 2 at large (3 in different QCD-based models JTj[ . 

Thus, we have established a region at large Q 2 and /3 ~ 1 where the diffractive 
structure function F®(/3, £, Q 2 ) is perturbatively calculable on the basis of the skewed 
gluon distribution. 



3 at realistic values of Q 2 and (3 

In this section, the diffractive cross section at high (3 and at realistic values of £ and Q 2 
is calculated on the basis of Eq. ([I]), corrected for skewing and the non-zero real parts of 
the amplitudes. Figure [I] shows the longitudinal structure function obtained with three 
different parametrizations of the leading order gluon distribution. 

In our calculation, we use a charm mass m c = 1.5 GeV and a leading order running 
coupling defined by ALo,n /= 3 = 144 MeV (a s (Mz) = 0.118) and evaluated at a fixed 
scale fi 2 = Q 2 /8. This choice is justified by the variation of e 2 , which determines the 
hardness of the process, between and Q 2 /4 in the massless case. Furthermore, we 
use the slope 6 = 8 observed in the closely related hard process of elastic p meson 
electroproduction [^1[ . 
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Figure I} The longitudinal structure function F£ ((3, £, Q 2 ) at Q 2 = 75 GeV 2 and £ = 
0.003 calculated on the basis of the leading order GRV |I7|,|(J (grv981o), CTEQ |]|[2(J 
(cteq51) and MRST |T9|,|20[] (mrs981o, central gluon) gluon distribution functions. The 
charm threshold is clearly visible. 



In the longitudinal case, no light quark masses are introduced. The endpoints of the z 
integration are cut off by demanding e 2 > s 2 min = 1.25 GeV 2 , which is the lowest scale at 
which the MRST gluon distribution is defined. If the cutoff is raised to e 2 min = 2.5 GeV 2 , 
the result changes by <5%. This is much less than other theoretical errors, e.g. the 
considerable uncertainty of the small- x gluon distribution (cf. Fig. |l|). 

In the following, we use the CTEQ gluon distribution because of its central position 
in the comparative plot of Fig. |l|. For this gluon distribution, the effects of skewing and 
of the real part at /3 = 1 are given by R 2 = 1.88 and 1 + C 2 = 1.42. 

Figure |2] shows our leading order results together with F 2 D data from HI. Unfortu- 
nately, at the values of f3 and Q 2 available at present, the transverse contribution cannot 
be neglected. Therefore we have also included an estimate of the transverse contribution 
in the plot. This contribution was calculated on the basis of Eq. (|l|) with skewing and 
real-part corrections. Unlike the prediction for F£ ', it is necessary to impose a regula- 
tor to estimate F® ■ Here we make the z integration finite by giving the light quarks a 
non-zero mass m q ~ A, which is justified by the production of massive hadrons in the 
final state. According to Eq. (f|), this mass gives rise to a cutoff z min and therefore to 
a minimal virtuality s 2 min = m 2 /(l — j3). At sufficiently large (3, this minimal virtuality 
is in the perturbative domain and a calculation on the basis of parametrizations of the 
gluon distribution becomes possible. We have plotted the corresponding curves in the j3 
range where e 2 min > 1.25 GeV 2 . 

Even though we have argued that e 2 is always in the perturbative domain, our 
transverse cross section should only be interpreted as an estimate. The reason for this 
is the unknown value of the regulator m q or, more generally, the lack of a quantitative 
understanding of the way in which confinement effects in the final state limit the possible 
kinematics of the qq pair. 
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Figure E|: The diffractive structure functions at Q 2 = 75 GeV 2 and £ = 0.003 (based on 
the CTEQ leading order gluon) compared with F<P data from HI f22fl . The estimates of 
the transverse contribution are based on mass regulators m q = 0.3 GeV (upper curves) 
and m q = 0.6 GeV (lower curves). 
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Figure [3|: Same as Fig. 0, but for Q 2 = 60 GeV 2 and £ = 0.0042 with corresponding 



data from ZEUS p2| 



In Fig. H|, the analysis of Fig. ^| is repeated for data from ZEUS. In both figures 
we have used the data extrapolated by the experimental collaborations to fixed £ and 
Q 2 values, because we are interested in the (3 shape at given £ and Q 2 . With more data 
points and smaller errors one could, in principle, extrapolate the data to (3 ~ 1 and thus 
directly compare with the perturbative calculation. 

At present, we can only say that the data are approaching the region where the 
calculable longitudinal contribution dominates. Our calculations, which include a skewing 
enhancement of almost a factor 2, are compatible with the data. Thus, the prospects for 
a quantitative description of _F 2 D data in this domain and for a model-independent test 
of skewing effects are good. 
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4 Beyond leading order 



At present there exists no complete next-to-leading order (NLO) calculation of diffrac- 
tive scattering. However, it is possible to improve the predictions beyond the leading 
approximation. Above we have already discussed the important contributions from the 
real part of the amplitudes and the effect from skewed parton distributions. Both are 
corrections to the leading log(l/x) approximation and lead to a significant enhancement 
of ■ Therefore, unless explicitly indicated, they will be taken into account (on ampli- 
tude level, see for details) in all the numerical results discussed below. In addition, 
the simple form of Eq. (|l|) was obtained in the leading \og{Q 2 /kj) approximation, where 
is the transverse momentum of the exchanged gluons that mediate the diffractive 
scattering. Only in this limit can the (loop-) integral over k? be performed immediately, 
leading to the integrated gluon distribution. In the more general case, the amplitudes 
contain a non-trivial /^-integration over the unintegrated gluon distribution /(£,&fOf] 



"Q 2 dfcf. 

This approach is an application of the so-called kr or high-energy factorization [24"| which 



£G(£,Q 2 ) = / rj: (13) 



is an extension of ordinary collinear factorization. It was already used in |23] for the 
prediction of diffractive open charm production and found to increase the cross section 
considerably. However, this effect is nearly compensated by the fact that the NLO gluon 
parametrizations, which are used for consistency in our NLO calculations, are smaller 
than the leading order ones. 

Two additional comments are in order here: when aiming at a full NLO prediction, 
in principle also qqg contributions have to be taken into account. While this is certainly 
true for arbitrary kinematics, additional real gluon emission is kinematically suppressed 
for large (3, which is the case we are interested in. Virtual gluon exchange, on the other 
hand, might be important. Up to now, such corrections to the ggqq vertex have not been 
calculated yet] 2 ] They might lead to an additional K factor of numerical relevance (see 
also |p3|1 ). 



Let us now discuss the numerical results of our NLO predictions for F® obtained 
in the /cy-factorization approach. Figure 4 shows the contributions to from F® and 
F® as continuous and dashed lines, respectively. Here Q 2 = 60 GeV 2 , £ = 0.0042, 
and we have used the CTEQ5M gluon as input in our calculations. The three dotted 
curves correspond to three different gluon parametrizations, as indicated in the plot; they 
demonstrate the strong dependence of our results on the gluon distribution used as input. 
For comparison we have shown two ZEUS data points, which are well described by the 
CTEQ5M prediction. We should, however, keep it in mind that our calculation is plagued 
by an uncertainty from the IR regime. The size of this uncertainty can be estimated by 
varying the light quark mass, which serves as an IR cutoff. Figure 5 shows predictions 



^or the full amplitudes including the integration over the gluon transverse momentum we refer the 
reader to 

2 In a first step the fmiteness of the impact factors was demonstrated, see |25f|. 
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Figure |4|: F 2 D (at Q 2 = 60 GeV 2 and £ = 0.0042) in NLO approximation as discussed in 
the text. The continuous (dashed) line shows the longitudinal (transverse) contribution 
calculated with the CTEQ5M gluon 



the dotted lines display the sum F® + F; 



D 
T 



for the three different NLO gluon parametrizations GRV [17| (grv98nlm), CTEQ [18|] 
(cteq5m) and MRST [Tj| (mrs99, central gluon) as indicated in the plot. ZEUS data |22 
are shown for comparison. 



for F® (using the CTEQ5M gluon) for the three different choices m q = 0/0.3/0.6 GeV.[] 
(In Fig. 4 m q = 0.3 GeV was taken as a default value.) Note that a cutoff given by the 
requirement e 2 + k 2 > (0.3 GeV) 2 , where K t is the transverse momentum of the produced 
quarks, together with m q = leads to results only slightly smaller than the upper curve 
of Fig. 5, obtained with m q = and no other 'confinement cutoff'. 

As long as (trivial) phase space suppression is not dominant, the variation of m q af- 
fects mainly the overall normalization. Together with the unknown K factor mentioned 
above this uncertainty at present makes it difficult to 'extract' the skewed gluon dis- 
tribution with high accuracy. Therefore, precise data at higher Q 2 and the calculation 
of the K factor are highly desirable. Nevertheless the strong sensitivity of F® on skew- 
ing effects is clearly demonstrated in Fig. 5 by the dash-dotted curve, which does not 
include skewing and has to be compared with the continuous line. Also note that the 
enhancement due to the real part contributions amounts to roughly 40 per cent, in good 
agreement with the enhancement factor found in the leading order analysis. 



3 Here we can calculate even for a vanishing quark mass. This is because no IR divergences appear 
in the integrals over the unintegrated gluon distribution, which is continued linearly for very small scales 
k\ < 1.5 GeV 2 as described in ||. 
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Figure §: F 2 D (at Q 2 = 60 GeV 2 and f = 0.0042) for three different choices of the 



light quark mass as indicated in the plot. For comparison two ZEUS data |22[] points are 
shown. The dash-dotted line is the prediction if the effects from skewing are neglected. 



5 Conclusions 



In the present letter, we have suggested a new way of extracting the small- a; skewed 
gluon distribution from diffractive electroproduction data. We have demonstrated that 
inclusive diffraction in the kinematic domain Q 2 3> M 2 ^> A 2 is, in principle, a per- 
turbatively calculable quantity, which is highly sensitive to skewing effects. Our leading 
order numerical analysis, which includes a skewing factor <2, is consistent with data 
now available in the region of interest. Next-to- leading order effects have been estimated 
using the unintegrated gluon distribution. The results reinforce the leading order situa- 
tion since the enhancement that goes with the explicit integration over transverse gluon 
momenta is largely compensated by the smaller next-to-leading order gluon distribution. 
It is clear that future, more precise data have the potential to constrain the small-x 
skewed gluon distribution significantly. 

Our method has the advantage that, unlike the case of diffractive vector meson 
electroproduction, no effects of non-perturbative final state wave functions affect the 
theoretical prediction. Furthermore, our calculation of inclusive diffraction at (3 — > 1 can 
be used to supplement the familiar leading twist partonic analyses of with a higher 
twist correction, which is known to be numerically important. Thus, we expect that the 
specific kinematic domain of F® at large Q 2 and (3 — > 1 will play an important role in 
the future development of hard diffraction, both theoretically and experimentally. 
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